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Abstract 



We investigate hadron production and attenuation of hadrons with strange and charm quarks 
(or antiquarks) as well as high transverse momentum hadrons in relativistic nucleus-nucleus col- 
lisions from 2 A-GeV to 21.3 A-TeV within two independent transport approaches (UrQMD and 
HSD). Both transport models are based on quark, diquark, string and hadronic degrees of freedom, 
but do not include any explicit phase transition to a quark-gluon plasma. From our dynamical 
calculations we find that both models do not describe the maximum in the /tt^ ratio at 20 - 
30 A-GeV in central Au+Au collisions found experimentally, though the excitation functions of 
strange mesons are reproduced well in HSD and UrQMD. Furthermore, the transport calculations 
show that the charmonium recreation hy D + D J/^ + meson reactions is comparable to 
the dissociation by 'comoving' mesons at RHIC energies contrary to SPS energies. This leads to 
the final result that the total J/^ suppression as a function of centrality at RHIC should be less 
than the suppression seen at SPS energies where the 'comover' dissociation is substantial and the 
backward channels play no role. Furthermore, our transport calculations - in comparison to exper- 
imental data on transverse momentum spectra from pp, d+Au and Au+Au reactions ~ show that 
pre-hadronic effects are responsible for both the hardening of the hadron spectra for low transverse 
momenta (Cronin effect) as well as the suppression of high px hadrons. The mutual interactions 
of formed hadrons are found to be negligible in central Au+Au collisions at y/s = 200 GeV for 

^ 6 GeV/c and the sizeable suppression seen experimentally is attributed to a large extent to 
the interactions of 'leading' pre-hadrons with the dense environment. 



1 Introduction 

The dynamics of ultra-relativistic nucleus-nucleus collisions at AGS, SPS and RHIC energies are of 
fundamental interest with respect to the properties of hadronic/partonic systems at high energy densities 
as encountered in the early phase of the 'big bang'. Especially the formation of a quark-gluon plasma 
(QGP) and its transition to interacting hadronic matter has motivated a large community for about 20 
to 30 years by now [1 . However, even after more than a decade of experiments at the SPS and in the 
last years at RHIC the complexity of the dynamics has not been unraveled and no conclusive evidence 
has been obtained for the formation of the QGP and/or the properties of the phase transition. 

* Supported by DFG 
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Apart from the light and strange flavor [u, u, d, d, s, s) quark physics and their hadronic bound states 
in the vacuum (vr, K, (p etc.) the interest in hadronic states with charm flavors (c, c) has been rising 
additionally in line with the development of new experimental facilities. This relates to the charm 
production cross section in pN, irN, pA and AA reactions as well as to their interactions with baryons 
and mesons which determine their properties (spectral functions) in the hadronic medium. Especially 
the J/^P suppression at SPS energies has been discussed as a signature for a QGP phase due to Debye 
screening pj. On the other hand, it has been pointed out - within statistical models - that at RHIC 
energies the charmonium formation from open charm + anticharm mesons might become essential [SI E] 
and even exceed the yield from primary A^A^ collisions. Thus charmonium suppression or enhancement 
is of central importance for an understanding of nucleus- nucleus collisions at RHIC energies. 

Additionally, transverse mass (or momentum) spectra of hadrons are in the center of interest. On 
the one hand the suppression of high transverse momentum hadrons is investigated in Au+Au reactions 
relative to pp colhsions at RHIC energies of a/s = 200 GeV 0, since the propagation of a fast quark 
through a hot colored medium (QGP) should be different than in ordinary high density hadron matter 
[HllZj. In fact, the PHENIX |H] and STAR [Hj collaborations have reported a large relative suppression of 
hadron spectra for transverse momenta pr above ~ 3-4 GeV/c which might point towards the creation 
of a QGP, since this suppression is not observed in d+Au interactions at the same bombarding energy 
per nucleon [TUIIIII- But it is not clear presently, to which extent this suppression might be due to 
ordinary hadronic final state interactions ^2); too- On the other hand, the measured transverse mass 
{rriT—) spectra of kaons at AGS and SPS energies show the opposite behavior for soft hadrons, i.e. a 
substantial hardening of the spectra in central Au-|-Au collisions relative to pp interactions (cf. also 

) . The flat spectra observed are commonly attributed to strong collective flow that is absent in the 
respective pp or pA data. An open question currently is, if the collective flow seen experimentally might 
be attributed to final state interactions of formed hadrons, too. 

In our studies we use two independent transport models that employ hadronic and string degrees 
of freedom, i.e. UrQMD ^1^] and HSD ^SlEI- They take into account the formation and multiple 
rescattering of hadrons and thus dynamically describe the generation of pressure in the early phase - 
dominated by strings - and the hadronic expansion phase. Both transport approaches are matched to 
reproduce the nucleon-nucleon, meson-nucleon and meson-meson cross section data in a wide kinematic 
range. We point out explicitly, that no parton-parton scattering processes are included in the studies 
below contrary to the multi-phase transport model (AMPT) 18 , which is currently employed from 
upper SPS to RHIC energies. Whereas the underlying concepts of UrQMD and HSD are very similar, 
the actual realizations differ considerably. 



2 Strange and charm hadrons 

2.1 Numerical implementation of 'hard' probes 

Whereas the dynamics of hadrons with u, d, s (and u, d, s) quarks is described nonperturbatively in 
HSD and UrQMD, the dynamics of 'rare' species needs a perturbative modeling in order to achieve 
acceptable statistics. For the production and propagation of open charm mesons, charmonia or hadrons 
with high transverse momentum (> 2.0 GeV/c) we employ the perturbative scheme as used in Refs. 
P^I^I^I22j . We recall that initial hard processes (such as cc, Drell-Yan or high pt hadron production 
from A^A^ collisions) are 'precalculated' to achieve a scaling of the inclusive cross section with the number 
of projectile and target nucleons as Ap x At when integrating over impact parameter (cf. Ref. [20]). 
Each open charm meson, charmonium or high pt hadron is produced in the transport calculation with 
a weight Wi given by the ratio of the actual production cross section divided by the inelastic nucleon- 
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Figure 1: The multiplicities 
for Ti+,r],K+,K-, 0, D,D and 
J/\E'-mesons for central colli- 
sions of Au + Au as a function of 
bombarding energy. In the HSD 
calculations elastic and inelastic 
reactions of all hadrons are in- 
cluded, but no in-medium modi- 
fications of their spectral func- 
tions (The figure is taken from 
Ref. EOI). 



nucleon cross section, i.e. 



= . (1) 

We then follow the motion of these hadrons within the full background of strings/hadrons by propagating 
them as free particles, i.e. neglecting in-medium potentials, but compute their collisional history with 
baryons and mesons or quarks and diquarks. For reactions with diquarks we use the corresponding 
reaction cross section with baryons multiplied by a factor of 2/3. For collisions with quarks (antiquarks) 
we adopt half of the cross section for collisions with mesons and for the leading quark (formed) baryon 
collision a factor of 1/3 is assumed. The final states in an inelastic collision are then modeled by the 
JETSET fragmentation scheme in the same way as for ordinary hadron-hadron collisions. 



2.2 Excitation functions of mesons in central collisions of heavy nuclei 

In order to provide an overview on meson production we show in Fig. 1 the calculated excitation 
function of 7c'^,K^,r],(f),D,D and J/\E' mesons in central Au + Au collisions from lower SIS to top 
RHIC energies without employing any self energies for these mesons in the HSD transport approach 
(cf. Ref. 120] )• We mention that the excitation functions of the strange hadrons from HSD and 
UrQMD match rather well with the data available, however, the sharp maximum in the /7r~^ ratio 
at ~ 20 to 30 A-GeV in central Au+Au (Pb+Pb) collisions is missed |23j. This also holds for the slope 
of the transverse mass spectra above ~ 5 A-GeV in central Au-|-Au (Pb-|-Pb) collisions |24[ I25j . 

The i)-mesons with a c are produced more frequently at low energies due to the associated production 
with Ac, Sc, similar to the kaon case, where kaon+hyperon production is more frequent than K + K 
production. At roughly 15 A-GeV the cross sections for open charm and charmonia are similar, while the 
ratio of open charm to charmonium bound states increases rapidly with energy. Again, this behaviour 
is comparable to the excitation functions in the strangeness sector when comparing K~^,K~ and 0- 
mesons. Since the excitation function for open charm drops very fast with decreasing bombarding 
energy, experiments at 20 to 30 A- GeV at the future GSI facility will be a challenging task since the 
multiplicity of the other mesons is higher by orders of magnitude. On the other hand, the perspectives 
for open charm measurements at RHIC appear promising since about 15 cc (or DD) pairs should be 
created in central Au + Au collisions at y/s = 200 GeV according to the HSD transport calculations. 
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2.3 Open charm and charmonia at SPS and RHIC 

We calculate open charm and charmonium production at SPS and RHIC energies (within the HSD trans- 
port approach for the overall reaction dynamics) using parametrizations for the elementary production 
channels including the charmed hadrons D, D, D* , D*, Ds, Dg, D*, D*, J/\E', \&(2S'), X2c from A^A^ and 
ttN collisions. The latter parametrizations are fitted to PYTHIA calculations [26j above y/s = 10 GeV 
and extrapolated to the individual thresholds, while the absolute strength of the cross sections is fixed 
by the experimental data as described in Ref. [201 • We here report on an extension of previous works 
|19t 123 I2Z| and include explicitly the backward channels 'charm + anticharm meson charmonia + 
meson' employing detailed balance in a more schematic interaction model with a single parameter or 
matrix element |Mp, that is fixed by the J/\E' suppression data from the NA50 collaboration at SPS 
energies (cf. Ref. p Bj ). 

Since the meson-meson dissociation and backward reactions typically occur with low relative mo- 
menta ('comovers') it is legitimate to write the cross section for the process mi + m2 + m4 
as 



where Ei and Si denote the energy and spin of hadron i, respectively, while Pi and Pf are the initial 
and final momenta for fixed invariant energy a/s. In Q \Mf\'^ stands for the effective matrix element 
squared, which for the different 2-body channels is taken of the form 



involving a single parameter Mq (cf. ||28j and Refs. therein). The relative factors of 3 in Q are guided by 
the sum rule studies in jSH] which suggest that the cross section is increased whenever a vector meson D* 
or D* appears in the final channel while another factor of 1/3 is introduced for each s or s quark involved. 
The factor ((M3 + M^)/^) in (0) accounts for the suppression of binary channels with increasing 
and has been fitted to the experimental data for the reactions tt + A^— >-p-|-A^, a;-|-A^, $ + A^, + A 
in Ref. [20] • The advantage of the model introduced in (j2)) is that detailed balance for the binary 
reactions can be employed strictly for each individual channel and the role of the backward reactions 
( J/\E'+meson formation hy D + D flavor exchange) can be explored without introducing any additional 
parameter once Mq is fixed. 

We directly step on with the results for the charmonium suppression and start with the system 
Ph + Pb at 160 A GeV to demonstrate that the 'late' comover dissociation model (j21) is approximately 
in line with the data of the NA50 Collaboration. The corresponding J/\E' suppression (in terms of 
the fi~^^~ decay branch relative to the Drell-Yan background from 2.9 - 4.5 GeV invariant mass) as a 
function of the transverse energy in Pb + Pb collisions at 160 A GeV is shown in Fig. 2. The solid 
line stands for the HSD result within the comover absorption scenario for the cross sections defined 
by (j2)) while the various data points reflect the preliminary NA50 data from the year 2000 (analysis 
A,B,C). A comparable agreement with the data is also achieved within the UrQMD model [23 (dashed 
histogram in Fig. 12)). We mention that there might be alternative explanations for J/\E' suppression 
as discussed in Refs. jnilE21 and/or further dissociation mechanism not considered here. However, for 
the purposes of the present study it is sufficient to point out that the cross sections employed here (cf. 
Figs. 6 and 7 in pHj) most likely are upper limits and do not lead to a sizeable recreation of charmonia 
by D + D channels at SPS energies. 
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Figure 2: The Jj^ suppres- 
sion (in terms of the de- 
cay branch relative to the Drell- 
Yan background from 2.9 - 4.5 
GeV invariant mass) as a func- 
tion of the transverse energy Ej' 
in Ph + Ph colhsions at 160 
A GeV. The sohd hne shows the 
HSD result within the comover 
absorption scenario presented in 
Section 2.3 |2H]. The different 
symbols stand for the NA50 data 
[3ni from the year 2000 (ana- 
lysis A,B,C) while the dashed 
histogram is the UrQMD result 
from Ref. 



For central Au + Au collisions at y/s = 200 GeV, however, the multiplicity of open charm pairs 
should be about 2 orders of magnitude larger than at 160 A-GeV (according to Fig. 1), such that a 
much higher J/\E' reformation rate (~ iV^g) is expected at RHIC energies (cf. Ref. In Fig- El we 

display the total J/\E' comover absorption rate (solid histogram) in comparison to the J/"^ reformation 
rate (dashed histogram) as a function of time in the center-of-mass frame for central Au+Au collisions 
at \/i = 200 GeV. The two rates become comparable for t > 4-5 fm/c and suggest that at the full 
RHIC energy of ^/s = 200 GeV the J/'^ comover dissociation is no longer important since the charmonia 
dissociated in this channel are approximately recreated in the backward channels. Accordingly, the J/^ 
dissociation at RHIC should be less pronounced than at SPS energies. Moreover, there is even a small 
excess of J/\E' formation hj D + D reactions in the first 2 fm/c. 

The preliminary data of the PHENIX Collaboration allow for a first glance at the situation 
encountered in Au + Au collisions at ^/s = 200 GeV. In order to compare with the preliminary data 
we have performed a rapidity cut Ay < 2 in the calculations. In Fig. |3 the J/\E' multiplicity per 
binary collision (times the branching ratio B) is shown as a function of the number of participating 
nucleons in comparison to the data at midrapidity. Since the statistics in quite limited so far on the 
experimental side, no final conclusion can presently be drawn, however, the data neither suggest a 
dramatic enhancement of J/"^ production nor a complete 'melting' of the charmonia in the QGP phase. 



3 High pt hadron suppression 

Before coming to actual calculations for deuteron-nucleus or nucleus-nucleus collisions it is important 
to briefly explain the concept of 'leading' and 'secondary' hadrons in the transport approaches, since 
this separation is of central importance for the results to be discussed below. To this aim we recall 
that in a high energy nucleon-nucleon collision two (or more) color-neutral strings are assumed to be 
formed, which phenomenologically describe the low energy 'coherent' gluon dynamics by means of a 
color electric field, which is stretched between the colored ends of each single string. The latter string 
ends are defined by the space-time coordinates of the constituents, i.e. a diquark and quark for a 
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Figure 3: The calculated rate 
of J/\l/ dissociation reactions 
with mesons (solid histogram) 
for central Au + Au collisions at 
y/s = 200 GeV in comparison to 
the rate of backward reactions of 
open charm pairs to J/\l/ + me- 
son (dashed histogram) accord- 
ing to the model specified in Sec- 
tion 2.3 (taken from Ref. |2H1)- 
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Figure 4: The calculated J/"^ 
multiplicity per binary collision 
- multiplied by the branching 
to dileptons - as a function of 
the number of participating nu- 
cleons Npart in comparison to 
the preliminary data from the 
PHENIX Collaboration ^ for 
Au + Au and pp reactions (taken 
from Ref. |2H|). 
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'baryonic' string or quark and antiquark for a 'mesonic' string. These constituents quarks (diquarks or 
antiquarks) are denoted as 'leading' quarks that constitute the 'leading' pre-hadrons. 

The time that is needed for the fragmentation of the strings and for the hadronization of the 
fragments we denote as formation time Tf. For simplicity we assume (in HSD), that the formation time 
is a constant Tf in the rest frame of each hadron and that it does not depend on the particle species. We 
recall, that due to time dilatation the formation time t/ in any reference frame is then proportional to 
the Lorentz 7-factor, i.e. tf = '-f ■ Tf. The size of Tf can be estimated by the time that the constituents 
of the hadrons (with velocity c) need to travel a distance of a typical hadronic radius (0.5-0.8 fm). 

We assume that hadrons, whose constituent quarks and antiquarks are created from the vacuum in 
the string fragmentation, do not interact with the surrounding nuclear medium within their formation 
time tf. For the leading pre-hadrons, i.e. those involving quarks (antiquarks) from the struck nucleons, 
we adopt a reduced effective cross section aiead during the formation time tf and the full hadronic 
cross section later on. Due to time dilatation light particles emerging from the middle of the string 
might escape the hadronic fireball without further interaction if they carry a high momentum relative 
to the rest frame of the fireball. However, the hadrons with transverse momenta larger than ~6 GeV/c 
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Figure 5: The ratio of 'leading 
particles' to 'all produced parti- 
cles' in A^ + A^ collisions {^/s=200 
GeV) at midrapidity for different 
particle classes function of 

transverse momentum within the 



PYTHIA model (taken from Ref. 
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predominantly stem from the string ends jHS] and therefore can interact directly with a reduced cross 
section (see below). 

The relative quark counting factors of 2/3 or 1/3 for the interactions of a 'leading' pre-hadron with 
a baryon or meson mentioned above might appear arbitrary and simplistic. However, this concept has 
been proven to work rather well for nucleus-nucleus collisions from SPS to RHIC energies 12H| as 
well as in hadron formation and attenuation in deep inelastic lepton scattering off nuclei [32] • Especially 
the latter reactions are important to understand the attenuation of hadrons with high (longitudinal) 
momentum in ordinary cold nuclear matter. The studies in Ref. |3S1 have demonstrated that the 
dominant final state interactions (FSI) of the hadrons with maximum momentum - as measured by 
the HERMES Collaboration - are compatible with the concepts described above. This also holds 
for antiproton production and attenuation in proton- nucleus collisions at AGS energies [3H]. Both 
independent studies point towards a hadron formation time Tf in the order of 0.4 to 0.8 fm/c. 

As mentioned above, the question of relevance is the fraction of leading pre-hadrons as a function of 
transverse momentum px for the different hadron species. This information is extracted from PYTHIA 
calculations and displayed in Fig. 5 for pp collisions at ^/s = 200 GeV. Here one notices slight differences 
between pions and especially antibaryons, however, the fraction of leading pre-hadrons increases with 
Pt and saturates above ~ 6 GeV/c. Thus at high momenta the major fraction of 'hadrons' is of leading 
quark, antiquark or diquark origin and - according to the assumptions stated above - may interact 
without delay with the partonic/hadronic environment. 

3.1 d+Au collisions at RHIC 

As pointed out before (cf. Ref. j21]) the slope of the transverse mass spectra especially of kaons 
and antikaons at SPS and RHIC energies is severely underestimated by the conventional transport ap- 
proaches HSD and UrQMD when including only hadronic reinteractions of secondary mesons. However, 
as known from the experimental studies of Refs. [321 HD] an enhancement of the transverse momentum 
cross section from p + A collisions relative to scaled pp collisions is already observed at ISR energies. 
This 'Cronin effect' is presently not fully understood in its details, but most likely related to an increase 
of the average transverse momentum squared < > of the partons in the nuclear medium prior to the 
'hard' scattering vertex. We speculate that such an enhancement of < > might be related to initial 
semi-hard gluon radiation of the leading parton in the medium, which is reduced in the vacuum due 
to color neutrality in a small volume. Since the microscopic mechanisms are beyond the scope of our 
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Figure 6: The suppression factor 
RdA (5) for minimum bias d + 
Au collisions {y/s=200 GeV) at 
midrapidity. Experimental data 
are from PHENIX (open 
circles) and STAR [H] (filled 
squares), while the hatched band 
shows the calculated charged 
hadron ratio (taken from Ref. 

)■ 



present analysis, we simulate this effect in the transport approach by increasing the average < > in 
the string fragmentation function with the number of previous collisions Nprev as 

< fc2 ^2 ^ aNprev)- (4) 

The parameter a ~ 0.25 — 0.4 is fixed in comparison to the experimental data for d+Au collisions 
[inillll (cf. Fig. 6). 

A comparison of the calculated ratio 

P f l/NZ'-'d'N,A/dydpT 
'^^^^^ < N^ou > /a;-/-^d^app/dydpT ^ ' 

is shown in Fig. IHlwith the respective data for charged hadrons from Refs. fOl ^] for a = 0.25 — 0.4 
(hatched band). In (0) < Ncou >~ 8.5 denotes the number of inelastic nucleon-nucleon collisions per 
event (for minimum bias collisions), whereas 0"^^^''^^ denotes the inelastic pp cross section 42 mb). We 
find that the HSD calculations give a rise in the ratio (0) for small px and a decrease for px > 2.5 GeV/c. 
The hatched band gives an estimate for the systematic uncertainty in Q with respect to the strength 
of the Cronin effect that carries over to studies for Au+Au reactions (see below). The important issue 
in this context is that no dramatic absorption of high pj- hadrons is found in the calculations as well as 
in the data. 



3.2 Au+Au collisions at RHIC 

We step on with 5% central Au+Au collisions at y/s = 200 GeV. For Au+Au collisions we define the 
nuclear modification factor in accordance with eq. (jSj) as 

P ( _ l/N^/X'"' d^NAA/dydpT 
^^^^"^^ (A^coii) /<i^^ d^apJdydpT ■ ^ ^ 

Fig. 13 shows a comparison of the calculations for eq. (jH)) with the data for charged hadrons from Refs. 
IHlin]- The surprising result in Fig. [3 is that the ratio Raa is roughly described in shape for the default 
HSD parameters (formation time Tp = 0.8 fm/c) within the uncertainties of the 'Cronin' simulations. 
Only the absolute magnitude of the high pt suppression seen experimentally is slightly underestimated. 
In order to understand this result we have decomposed the ratio Raa into contributions from secondary 
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Figure 7: The suppression fac- 
tor Raa (6) of charged hadrons 
at 5% (or 10%) central Au + 
Au colhsions (v/i=200 GeV) at 
midrapidity for a = 0.25 to 
0.4 (hatched band). The ex- 
perimental data are from Refs. 
jHl Ej and show some additional 
suppression mechanism, which 
might be attributed to partonic 
interactions in a colored medium 
(taken from Ref. [35J). 



pions, kaons and p,uj,K* vector mesons (before decays) as well as from leading pre-hadrons. We find 
that formation time effects play a substantial role since heavier hadrons are formed earlier than light 
pions in the cms frame at fixed transverse momentum due to the lower Lorentz boost. However, the 
total suppression by scatterings of formed hadrons is negligible for Pt ^ G GeV/c and only of relevance 
at lower pt for the shape of Raa{pt)- Furthermore, the attenuation of the leading pre-hadrons is 
roughly independent on px and hadron type and gives Raa ~ 0.1. Thus a large part of the attenuation 
seen experimentally should be addressed to the interactions of leading pre-hadrons with already formed 
hadrons. Nevertheless, the experimental data indicate some additional suppression mechanism. 

4 Conclusions 

Summarizing this contribution, we point out that strange hadron production in central Au+Au (or 
Pb+Pb) collisions is quite well described in the independent transport approaches HSD and UrQMD 
The exception are the pion rapidity spectra at the highest AGS energy and lower SPS energies, 
which are overestimated by both models. As a consequence the HSD and UrQMD transport approaches 
underestimate the experimental maximum of the /tt^ ratio at ~ 20-30 A-GeV |i231 125j . 

Furthermore, the transport calculations show that the charmonium recreation by D -|- Z) ^ J/^ + 
meson reactions is comparable to the dissociation by 'comoving' mesons at RHIC energies contrary to 
SPS energies. This leads to the final result that the total J/\l/ suppression as a function of centrality 
at RHIC should be less than the suppression seen at SPS energies where the 'comover' dissociation is 
substantial and the backward channels play no role. Since the statistics is quite limited in the present 
experiments for the J/\l/ yield versus centrality, no final conclusion can be drawn so far; however, the 
available data neither suggest a dramatic enhancement of J/\l/ production nor a complete 'melting' of 
the charmonia in a QGP phase [2H1- 

Moreover, our transport calculations in comparison to experimental data on high transverse momen- 
tum spectra from d+Au and Au+Au reactions show that pre-hadronic effects are responsible for both 
the hardening of the spectra for low transverse momenta J^T (Cronin effect) as well as the suppres- 
sion of high Pt hadrons [SSj • The interactions of formed hadrons are found to be negligible in central 
Au+Au collisions at ^/s = 200 GeV for > 6 GeV/c and the large suppression seen experimentally 
is attributed to a large extent to the interactions of 'leading' pre-hadrons with the dense environment, 
which should be partly of partonic nature in order to explain the large attenuation seen in central 
Au+Au collisions. 
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